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Abstract 

We present a simple proof of the all-order exponentiation of soft logarithmic corrections to 
hard processes in perturbative QCD. Our argument is based on proving that all large logs 
in the soft limit can be expressed in terms of a single dimensionful variable, and then using 
the renormalization group to resum them. Beyond the next-to-leading log level, our result is 
somewhat less predictive than previous all-order resummation formulae, but it does not rely on 
non-standard factorization, and it is thus possibly more general. We use our result to settle 
issues of convergence of the resummed series, we discuss scheme dependence at the resummed 
level, and we provide explicit resummed expressions in various factorization schemes. 
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1 The structure of soft resummation 



The resummation of large logarithms related to soft radiation near the boundary of phase space 
is a classic subject, which stems from very early results in quantum electrodynamics |l|], and is of 
great phenomenological and theoretical interest for QCD at colliders . The exponentiation of 
soft logs related to gluon emission has been proven in QCD to leading order 0, and to next-to- 
leading order (NLO) through an explicit analysis of the relevant Feynman graphs, either using 
eikonal or factorization techniques, which lead to equivalent ^ resummation formulae. 
These resummation formuale were subsequently used |0] to obtain resummed results for physical 
cross sections. 

More recently p|, it was shown that a large class of resummation formulae can be derived 
in a unified way and to all orders by assuming the validity of a suitable factorization ||^, and 
then applying the renormalization group to it. The price to pay for this generality is that 
one has to rely on a rather nontrivial two-scale factorization, which goes beyond the standard 
factorization of coUinear large logs. Also, the resummed result obtained in this way looks quite 
formal: in particular, its equivalence with resummation formulae of the form of Refs. @, ^ is 
only estabhshed to next-to-leading order. 

In this paper, we present a simple argument for the all-order resummation of soft logs. Our 
argument is based on the following steps. First, we establish to all orders the relation between the 
resummation of large logs of the kinematic variable 1 — x, and large logs of its Mellin conjugate 
variable N, up to power-suppressed terms in the soft limit x ^ 1, i.e. N ^ oo. This relation 
generalizes to all orders the relation between the next-to-leading resummations of Refs. |^ and 

which can be respectively viewed as a resummation of ln(l — x) contributions to the splitting 
function, and a resummation of In contributions to the anomalous dimension. 

Next, we prove, that the large logarithms of (1 — x) all stem from logs of a dimensionful 
variable — x)"', where a is an integer which characterizes the physical process. This origin 
of this result is essentially kinematic. It is related to the structure of the phase space for tree-level 
n-particle emission, which turns out not to be modified by loop corrections. Finally, we resum the 
logs of this latter dimensionful variable using standard renormalization-group techniques. This 
leads to a resummed result which generalizes to all orders the next-to-leading log resummation 
of Refs. 1^, H in a way which is somewhat less predictive than the all-order result of Ref. 0. 
Indeed, we also conclude, like Ref. [|], that logs of N exponentiate to all orders. However, using 
the resummation formula of Ref. [§] the N''~^LO resummed result is determined by a fixed-order 

N'^LO calculation, whereas using our resummation the N'^~^LO resummed result is determined 

fc(fc+i) 

by a fixed-order N 2 ^^LO computation. On the other hand, our result does not require 
two-scale factorization, and it is thus easier to derive, and possibly more general. 

We then exploit our relation between the resummation of logs of N and 1 — a; to study the 
relation between various resummed results. First, given that the resummed result can be equiv- 
alently viewed as a ln(l — x) resummation of the splitting function or a In resummation of the 
anomalous dimension, we discuss the possibility of further casting it as a ln(l — x) resummation 
of the physical cross section. We show that, whereas this is formally possible, the ensuing result 
is meaningless since it is the sum of a factorially divergent 0] series. We trace the origin of this 
divergence, and we show that it is related to subleading logs, and thus in particular unrelated 
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to power corrections. We then turn to the derivation of exphcit expressions for the resummed 
cross section in a variety of factorization schemes, and discuss the wider freedom in choice of 
factorization which is available at the resummed level. 

2 Next-to-leading log resummation 

In this section, we review the well-known resummation formulae of Refs. @, ^ , and recast them 
in a form which is best suited to our subsequent derivation and generalization. 

We consider the total cross section for a partonic hard process in the vicinity of the kinematic 
boundary for the production of the final state. We assume that the kinematics of the process 
is specified by a hard scale Q^, and a dimensionless variable x, defined in such a way that 
< X < 1, and x = 1 is the kinematic boundary, i.e. the limit in which there is no phase space 
for the emission of extra particles. Kinematically, the processes we consider fall into two broad 
categories, according to whether the large scale which makes the process perturbative is the 
virtuality of one of the incoming particles, as in deep- inelastic scattering (DIS), or the invariant 
mass of the final state, as in Drell-Yan (DY). We will not consider processes, such as jet or 
prompt-photon production, in which the hard scale is provided by a transverse energy. As is 
well known, the radiation of extra soft partons close to the kinematic boundary enhances the the 
total cross section by contributions proportional to powers of ln(l — x). These logs are generated 
by the integration over the phase space of the (on-shell) emitted partons, and therefore, since 
this phase space has an azimuthal symmetry, we expect an extra emission in the final state to 
be accompanied by at most two powers of the large log. 

Consider specifically the Mellin moments of the cross section 

a{N,Q^) = [^dxx^-^a{x,Q^) (2.1) 
Jo 

and assume that for DIS-like processes they can be written in factorized form as 

a(iV,Q2) = c(iV,QV/i',«,(/i2)) F(iV,/i2) 

= C(iV,l,«,(g2)) F{N,Q'), (2.2) 

where F{N, /i^) are Mellin moments of a parton density F{x, /i^), C{N, Q"^/ /i^, Q;s(/i^)) is a hard 
perturbative coefficient, and the last equation gives the standard renormalization-group improved 
form of the factorized cross section. In the case of Drell-Yan-like processes, there are two partons 
in the initial state, and thus Eq. ( |2.2| ) should be modified accordingly. In general, F in Eq. ( |2.2| ) 
is an array of parton distributions (quarks and gluons), but the ensuing complication of the 
formalism is trivial. Furthermore, to next-to-leading order, as we will see explicitly in the sequel, 
there is no quark-gluon mixing, so F can be effectively taken to be a single parton distribution. 
The validity of a factorization of the form Eq. ( |2.2D is the only assumption we will make on the 
physical processes to which our resummation applies, and it amounts to the assumption that 
the hard coefficient C {N,Q'^/fi'^,as{fJ''^)) can be multiplicatively renormalized. We will discuss 
this in more detail in Sect. |4.2| . 
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Because resummation takes the form of an exponentiation, it is convenient to classify large 
logs in terms of the log derivative of the cross section a, i.e., the so-called physical anomalous 
dimension ITOll defined as 



Q 



,da{N,Q^) 



^{N,as{Q'))a{N,Q'). 



(2.3) 



When addressing issues of operator mixing, it is convenient to choose a basis of physical cross 
sections (or physical observables) equal in number to the independent parton distributions, so 
the physical anomalous dimension is a matrix. The physical anomalous dimension 7 Eq. (|2.3| ) is 
independent of factorization scale, and it is related to the standard anomalous dimension 7"^^, 
defined by 

>dF{N,fi^) 



7^^(iV,a,(;.^))F(iV,;x^), 



according to 



(2.4) 



(2.5) 



It follows that the physical and standard anomalous dimensions coincide at leading order, but 
differ beyond leading order. In terms of the physical anomalous dimension, the cross section can 
be written as 



where 



a{N,Q') = KiN;QiQ')a{N,Q',) = exp E{N;QiQ') a(iV,Q^) 



(2.6) 



(2.7) 



J . ^7^^(iV,«,(P)) + lnC(iV,l,a,(Q2)) -lnC(iV,l,a,(Q^)). 



The structure of leading logs is then found by noting that upon Mellin transformation 
Eq. ( |2.1|) the large x region is mapped onto the large N region, so logs of 1 — x become logs of 
N and conversely. Explicitly 



dx X 



N-l 



- x) 
1 — X 



dx 



X 



N-l 



1 — X 



- x) 



p+l 



lnP+^ — - 7r In^ ^ + O InP 



N 



N 



N 



(2.8) 



where je is the Euler gamma, i.e., je = ""^(1)- The resummed next-to-leading order result of 
Refs. H, |5I1 has the form 



dx 



X 



N-l 



1 — X 



(2.9) 



where 

gi^s) = gias + g2a1, (2.10) 

and a = 1 for processes with deep-inelastic kinematics, and a = 2 for Drell-Yan kinematics. This 
result implies that the physical anomalous dimension is given by 

„1 ™Af-l _ 1 

'j{N,as{k^)) = a dx — \gias{k\l - x^) + g2a^{e{l - x)^) 



1 — X 

+0 (a,^+'(A;2) In^ A^) + 0(iV°) (2.11) 
= 7iK(A:') IniV) + a,(fc2)72(a.(A:') IniV) + O (^^'(A;^) In'^ n) + 0(iV°), 

where 0(A^°) denotes non- logarithmic terms. 

This means that, up to terms which are finite as ^ oo, the physical anomalous dimension 
can be expanded as a power series in at fixed In A^. The leading log contribution to the 
physical anomalous dimension is the sum of terms to all orders in as where the powers of as and 
In 1/N coincide; the next-to-leading order terms have one more power of as, and so on. Because 
the QCD beta function starts at order a^, the order a^ contribution to 7 Eq. ( p.ll|) leads to 



a contribution of order a^(Q^) In'^"'"^ to E Eq. (|2.7] ). In particular, the 0{as) contribution 
corresponds to a series of double- log contributions to the cross section (in InA^), in agreement 
with the simple power counting which leads to expect at most two extra powers of In A^ when 
the cross section is computed at an extra order in as- Notice, however, that the power counting 
for the physical anomalous dimension Eq. (|2.11|) is in fact more restrictive than what allowed 



on the basis of the fact that the cross section can contain at most double logs, in that it is free 
of contributions of order aj"*"^ In^'^ A^ with k > 1. Indeed, Eq. (|2.11|) implies that the anomalous 
dimension contains at most single logs, i.e., each extra power of as is accompanied at most by 
an extra power of In A^. 

The resummed result Eq. ( p.ll|) implies that all the leading-log coefficients are determined 



in terms of the single parameter gi, and the leading and next-to- leading ones by gi and g2. at 
leading log we have 

^{N,as{Q^)) = ag, / dx — «,(Q2) ^(_/3o)i «i(Q2) jj^.^l _ ^^^a^ ^2.12) 

Jo I - X 

where we have expanded the beta function as 

fi'^asifi') = f3{as) = -/3oa'(/^') - Pia^ifi^) + 0{at). (2.13) 



The term with j = in the series Eq. ( p.l2|) is a contribution to the leading-order anomalous 



dimension: hence, gi is just the coefficient of the ln(l/A^) contribution to the usual one-loop 
anomalous dimension. Similarly, it is apparent that ^2 is the coefficient of the In A^ contribution 
to the two-loop anomalous dimension. 

If Eqs. ( p^ , |2.1CI|) hold beyond next-to-leading order, as claimed in Ref. p, i.e., if g{as) is 



a power series in as with numerical coefficients gk, then gk is just the coefficient of the InA^ 
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contribution to the fc-loop physical anomalous dimension: the next'^-to-leading log resummation 
is determined by knowledge of fixed-order [k + l)-loop result. In Sect. ^ we will prove a some- 
what less predictive form of the resummation, where the next'^-to-leading log resummation is 
determined in general only by knowledge of a higher fixed-order computation. 

Beyond leading order the standard anomalous dimension differs from the physical one, so gk 
with k > 1 receives a contribution both from the standard anomalous dimension and from the 
coefficient function. It is thus natural to rewrite the resummation formula Eq. ( p.9|) in analogy 
to the unresummed separation Eq. (|2.7|) of contributions to E from the anomalous dimension 
and coefficient function, by separating off the contribution which originates from the anomalous 
dimension 7^^ Eq. (|2] 



1 x^-i - 1 



a dx 
Jo 1 



X 



Hi-xr dk"^ 



+ BMQ^l - xy)) - BMQlii - xr)) 



(2.14) 



where A{as) = Aias + ^2^^ is defined as the sum of one and two-loop coefficients of ln(l/A^) in 
the standard anomalous dimension, and 



B{as) = Bias 

makes up for the difference between physical and standard anomalous dimensions: 

dBiasik')) 



g{as) = A{as) + 



din fc2 



We can then rewrite the resummed cross section 

a'^^'iN, Q2) = exp Ie'^^^N; QI, Q^)] a'^^N, Ql) 



(2.15) 



(2.16) 



(2.17) 



in factorized form according to Eq. ( p.2| ) by collecting all Q^-dependent contributions to the 
resummation Eq. ( p.l4|) into a resummed perturbative coefficient C'^'^^ 



a-^(iV, Q2) = ^-(^^ g2/^2^ a,(/x2))F(iV, /x^) 



(2.18) 
(2.19) 



1 ^N-l _ I 

exp < a I dx 

^ Wo l-x 



^Aiasik')) + BiasiQ\l-xr)) 



The precise definition of the parton distribution F and the factorization scale /z^ will depend on 
the choice of factorization scheme: according to the choice of scheme, the resummed terms will 
be either part of the hard coefficient C^^^, or of the evolution of the structure function F. A 
detailed discussion of factorization scheme choices will be given in Sect. ^. 

The resummed results Eq. (|2.9| ), or Eqs. (|2.18[ - ^n9D , can be explicitly cast in the form of 
an exponentiation of logs of 1/A^. To this purpose, we note that Eq. (|2.8|) implies that, up to 
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next-to-leading log level 
Jo 



dx 



X 



1 — X 



lnP(l - x) 



1 



p+l 

"1- 



P+i 



(2.20) 







dx 
1 — X 



1 - 7s 



d 



d\n{l — x) 



lnP{l-x). 



The derivative term in the last expression is the first subleading correction to the leading-log 
term. Note that, whereas a next-to-leading log term can always be written as a derivative of 
the leading- log result, it is nontrivial that the coefficient of this derivative in Eq. ( p. 201) does not 
depend on the power p of ln(l — x). This fact can be exploited to rewrite the resummed result 
Eq. (|19D as 



lnC"^^^(iV,gV/i^a,(/i^)) 



dx 





dn 
1 n 



1 — X 



^Aia,ik'))+Bia,iQ'il-xr)) 



Q'dll 

2 f^2 



A{as{kyn)) + B{a,{Qyn)) 



(2.21) 



where in the second step we have performed the change of variable n = {1 — x) °, and to this 
order B{as) = B{as) — a'^E^iC^s- In Eq- ( |2.21| ) the exponentiation is directly expressed in terms 
of logs of 

Similarly, we can express the resummed physical anomalous dimension Eq. ( |2.9| ) directly in 
terms of logs of by using Eq. ( |2.20|) to perform the Mellin transform. We get 



1-^ dx fQ' dk^ 



1-x Jo?, A;2 



■gMk\l-xr)) 



dn fQ' dk'^ 
1 n JqI k^ 



where 



1 - ajE 



d 



dlnfi'^ 



g {as{k'^/n)^ 



g{as{n^)). 



(2.22) 
(2.23) 

(2.24) 



The physical anomalous dimension Eq. ( |2.3| ) which corresponds to the resummed result 
Eq. ( p. 231) is particularly simple: we get 



7(iV,a,(Q^)) 



dn 



1 n 



giOisiQ /n) +g2a^{Q /n 



g, asiQVN-) 1 ( ' 



(2.25) 
(2.26) 



where the coefficients of the two expansions Eq. (|2.11|) and Eq. ( p.25| ) up to next-to-leading order 
are related by Eq. ( p.24| ): 



91 



92 



(§2 - a'jEpogi) ■ 



(2.27) 
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The leading- and next-to-leading-log coefficients in the expansion of Eq. (|2.11 



) are given by 



7i = - 



72 



Po 



l + /3oa.(Q')ln — 



In 



l + /9oas(Q' 



In — 



l + /5oa.(Q2)ln^ 



(2.28) 
(2.29) 



where we have used the beta function Eq. ( p.l3|) . 

Summarizing, the resummation of Refs. 0, Q can be compactly expressed in the two equiva- 
lent forms Eq. ( p.llj ) or Eq. (|2.25| ) of the resummed physical anomalous dimension. The former 
shows that the resummed anomalous dimension is the Mellin transform of a resummed leading-, 
next-to-leading-,. . . ,ln(l — x) splitting function, which in turn is a power series in as(Q^(l —x)"") 
with numerical coefficients. The latter shows that this Mellin transform can be just expressed 
as a series in ^^(Q^/A^"), with numerical coefficients. Knowledge of the fc-th coefficient of either 
expansion determines the next*^" ^-to-leading resummation. This coefficient can be read off the 
coefficient of the In N term in a fixed fc-th order computation of the physical anomalous dimen- 
sion. Up to next-to-leading order, the two forms of the resummation are related by Eq. ( p. 24 ). 

Our aim here is to derive an all-order generalization of these resummation formulae. A 
preliminary step is thus to obtain an all-order generalization of the relation between the two 
available forms of the next-to-leading log resummation. This generalization will be derived in 
the next section. 



3 N space and x space beyond next-to-leading order 

We wish to find an all-order generalization of the relation between leading InA^ and leading 
ln(l — x) resummation discussed in the previous section. The leading InA^ resummation is 
expressed in terms of 

= \n^^ = -p lnP-^(l - x), (3.1) 



^ ^ Jo l-x 

while the leading ln(l — x) resummation is expressed in terms of 

J„ = / dx lnP(l - x). (3.2) 

Jo 1 — X 

In order to relate the two resummations we must therefore generalize to all orders the next- 
to- leading order relation Eq. ( |2.20[ ), i.e., compute the Mellin transform Eq. ( |3.2D to all-order 
logarithmic accuracy, up to power corrections. 

To this purpose, we notice that all integrals Ip can be obtained from a generating function 
Giv): 

G{7]) = r dx {x''-^ - 1) {1 - xy-\ (3.3) 
Jo 
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The integral can be determined at large using the Stirling formula: 



T{N)T{t]) 1 1 



r(i + r/) 



r(iv + ?7) 7] T] 

On the other hand, the generating function for Lp is just N~^: 

dP 



O 



(3.4) 



Lp = — iV-" 
^ dT]P 



(3.5) 



ri=0 



Hence, Eq. ( p.4|) can be viewed as a relation between the generating functions for Ip and Lp. 
In particular, Taylor-expanding r(l + rj) in the expression Eq. (p.4|) of G{ri) leads to leading, 
next-to-leading, ... In relations: 



1 oo 

— y 

k=0 



^ r('=)(i) rf^' 

to k\ din' Jo 



dx(l-x)'^-\ 



(3.6) 



For the individual Ip Eq. ( |3.3|) we get 
-'p ~ ~ 



fe=o k\ din" Jo 



1 — X 



N 



Because 



d' 



d In' ^ ^0 



dx 



InP(l-x) /-i-^ rfx d'ln^il-x) 



1 — X 



Eq. ( p.7|) can be rewritten as 



-E 

fc=0 



d' 



l-x dln^d - X 



l-x dln'{\-x) 



■lnP(l-x)+r(^'+^)(l)+0 



(3.7) 



(3.^ 



(3.9) 



This is the desired generalization of Eq. ( p. 201 ) to all orders. Notice that the coefficient in the 
expansion remain p-independent to all orders. Performing the derivatives explicitly, we get 



p+i 



V + 



fc=0 



p + 1 
k 



r<'.(i,(in- +o(-|. 



(3,10) 



The inverse result, expressing Lp in terms of /p, can be analogously found by inverting the 
relation between generating functions: 



riG[r]) + 1 

r(i + r7) ■ 



(3.11) 
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Proceeding as above, we get 

ln"± = yfl kl^ / dx- ^ in"(l-x) + A(")(l) + 0(-), (3.12) 



where A^''\r]) is the A;-th derivative of 



MV) - (3.13) 



Evaluating the integrals explicitly we get 



In" ^ = E ( ^ J A(-)(l) z + AW(1) + O (^) . (3.14) 



Because of the p-independence of the coefficients of the expansion Eq. (|3.9|) , we can determine 
explicitly the Mellin transform of a generic function 

oo 

gi\nil-x)) = J2gp\nPil-x) (3.15) 

p=0 

of ln(l — x), up to power corrections: 

/ dx— g {\n{l - x)) = J2 9p dx — InP(l-x) 

Jo 1 — X ^ JO 1 — X 

1 ^ -^^MHi - .)) + o,;v») (3.16) 

r-l fc! Jo 1 - a; rfln (1 - x) 



fc^o /c! Jo 1 - a; ciln^Cl - x 

'•'^-w dx 
lo 1 — 



/""-^^(ln(l-x))+0(iVO) 

JO 1 — X 



where 



,a„(i-.)).-g£^-^^-|-^j(Ki-.)). (3.17) 

and the last (constant) term in Eq. ( |3.9| ) has been dropped, so Eq. ( |3.16| ) is only correct up to 
non- logarithmic terms, denoted by 0{N^). 

The inverse relation can be analogously derived. Namely, we can cast the integral of any 
function 

oo 

g{ln{l-x)) = J2g,lnP{l-x) (3.18) 
as a Mellin transform, up to non-logarithmic terms: 



dx- gilxiil - x)) + OiN"^) (3.19) 

1 — X 
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where in the first step we have used Eq. (|3.12 



), and now 



^(ln(l - x)) 



g A(^)(l) 



k=0 



k\ dlifil- 



X] 



-g{\\i{l - x)). 



(3.20) 



Equations ( p.l6| - |3T20| ) immediately provide us with the sought-for all-order generalization 
of the relation between the two forms of the resummed anomalous dimension Eq. (|2.11|) and 
Eq. ( ^.251 ), or of the resummation factor, Eqs. ( f^.9D and ( |2.23| ). Specifically, 



1 , x^-i - 1 

dx 

1 



i-N'^ fin 

g[as{k\l-xT)]= / —g[a,{e/n)] 



with 



p=0 



r(p)(l)aP dP 
p\ d In^ /i^ 



(3.21) 



(3.22) 



By choosing a factorization scheme, we can then obtain from the resummation of the physical 
anomalous dimension a resummation of the coefficient function e.g. of the form of Eq. ( |2.19|) , 
with the functions A and B computed to the corresponding order. 

For future applications, it is interesting to observe that, through similar manipulations, it is 
also possible to derive a relation between a generic function of In and a function of ln(l — x). 
Indeed, given a function 



h In 



Eq. ( gl2D implies that 



h In 



N 



dx 



X 



N 



N-l 



J2h,lnP 



N 



(3.23) 



d^ 



1-x dln^il-x 



-h{\n{l - x)) + 0{N^ 



(3.24) 



The right-hand side of Eq. ( p. 241) can be viewed as the Mellin transform of a function (more 
properly a distribution) h(\n{l — x)): 



dx x^ ^h{ln(l — x)); 



/i In — 

V N 

oo /^(fe-l)fn 



d'' 



1 -xdln^d - X 



-f{\n{l-x)) 



+ 0(A^°). (3.25) 



4 All-order resummation 



This section contains the main result of our paper, namely the proof of an all-order generalization 
of the next-to- leading order resummations discussed in Sect. This proof exploits the relation 
between ln(l — x) resummation and InA^ resummation established in Sect. |^, and proceeds in 
two steps. First, (in Sect. [4.1|) we determine the N dependence of the regularized coefficient 
function in the large- iV limit. We show that the kinematic structure of the fc-particle phase 
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space essentially determines the form of the iV-dependence in the large limit, and that loop 
integrations do not modify this result. Then (in Sect. [4.2| ) we prove that, given this form of the 
iV-dependence of the regularized cross section, renormalization group invariance fixes the all- 
order dependence of the physical anomalous dimension Eq. ( p.3| ) in such a way that an infinite 
class of leading, next-to-leading,. . . resummations are found in terms of corresponding fixed order 
results. In particular, in Sect. 13 we show that the next-to-leading- log resummation Eq. (|2.23|) 
can be obtained by exploiting available fixed-order results. 



4.1 Kinematic structure of soft logs 

We consider the perturbative expansion of the bare coefficient function in powers of the bare 
coupling constant ao, 

oo 

CW(x, Q^ ao, e) = E (x, Q^ e) (4.1) 

n=0 

and its Mellin transform 

oo 

C(°)(iV,Q^ao,e) = ^ «o (A^, Q', e) (4.2) 

n=0 

in dimensional regularization with d = A-~2e space-time dimensions. By coefficient function here 
we mean a vector whose elements are the parton-level cross section for the partonic subprocesses 
that contribute to the given process: the formalism and the results of this section remain true 
even in the presence of parton mixing. 

The result that will be proven in this section, and then renormalization-group improved in 
the next section, is that in the large- limit, C^^^ (N, Q^, e) has the following structure: 

C(^\N,Q\e) = ±d:i{e){Qr^--'^' (|^) '\o{l/N), (4.3) 

where 0{1/N) denotes all terms which vanish as ^ cxo in the limit e — > 0. We will also see 
that the coefficients C^^ij have a pole of order 2n in e, related to infrared singularities. Even 



though such poles cancel in the coefficient function, which is free of infrared singularities |11 



the interference of the poles with the powers of A^~^ is responsible for the presence of powers 
of In N in the renormalized four- dimensional cross section. Therefore, soft logs may be viewed 
as being due to an incomplete cancellation between real and virtual contributions to C^^ in the 
soft limit. 

Equation ( [4.3|) will be established by proving that 



(x,g^e) = (g2)- 



n -^(0) 

Ci^^ie) 6(1 -x) + Y: (1 - x)-'-^'^ 



+ (l-x)° , (4.4) 



where O [(1 — x)°] denote terms which are not divergent as x — 1 in the e — limit. Eq. (O) 
implies Eq. (^4.3| ) because 

, ' dxx^^' (1 - x)-^-"'^^ = ^(^)^(~^^;^) = Y(-ake) N"''' + 0(1/N), (4.5) 
r(A^ — ake) 
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while the Melhii transform of any function of x which is not divergent as a; — > 1 {e.g. any 
non- negative power of 1 — a;) vanishes as — * oo. 

The content of Eq. (|4.3| ) is that, in the large N limit, the dependence of the regularized 
cross section on only goes through integer powers of the dimensionful variable (Q'^/N^) 
Because in d dimensions the coupling constant has the dimensions of Q^^, the dependence on 
this dimensionful variable is related through the renormalization group to the running of the 
coupling. As we shall show in Sect. [4.2| , this is sufficient to establish the desired all-order 
resummation. In the remainder of this subsection we will show that Eq. ([4.4|) essentially follows 



from the kinematics of soft emission. The reader who is not interested in technicalities can skip 
directly to Sect. where only the result Eq. (|4.3| ) will be used. 

We consider deep-inelastic scattering and Drell-Yan production as representative cases of the 
two classes of processes discussed in Sect. ||. In the case of deep-inelastic scattering, the relevant 
parton subpro cesses are 

7*(g) + Q(p) ^ Qip) + X (4.6) 

l*{q) + Q{p) - Q{p) + X, (4.7) 

where Q is a quark or an antiquark, Q a gluon, and X is any collection of quarks and gluons. 
In the case of Drell-Yan, the relevant subprocesses are 

Q{p) + Qip) ^ l*{q) + X (4.8) 
Q{p)+Q{p') ^ l*{q) + X (4.9) 
Q{p)+Q{p') ^ 7*{q) + X. (4.10) 

We are interested in the most singular contributions to the cross section as x — 1, and 
specifically in the way these terms depend on 1 — x. To this purpose, we consider the dependence 
on 1 — X of the phase space and of the amplitude for the generic processes Eq. (0- 4.10 ). The 
structure of the amplitude will in turn be discussed by considering first tree-level processes, and 
then processes with loops. We will first carry out the proof for deep-inelastic scattering, and 
then for Drell-Yan. 

As discussed in the Appendix, the phase space for a generic process with incoming momentum 
P and n bodies in the final state with outgoing momenta ki,. . . ,kn can be expressed in terms 
of two-body phase space integrals by using recursively the identity Eq. ([A.5|) . For DIS-like 
processes, all outgoing particles are massless partons. Indicating with p and p' the incoming and 
outgoing quark momenta and with q the 7* momentum we have 

r dMl 



+ q;ki,...,k 

niP ) 



27r 



dhip + q; kn, Pn) d(j)n{Pn] ki 



where 



{p + qf 



Q\i-x) 



X 



Applying Eq. (|4.11 

d4>n+l{Pn+l'-, kl, 



recursively we get 
• ; kn, p) = 



d(f)2{p + q; k, 



2tx 
Mi dMi 

2n 



ni Pn) 



Ml dM: 



n—l 



27r 



.XP3;k2,P2)dMP2;k^,P, 



kn-up'), (4.11) 
(4.12) 

2{Pn', kn-1, Pn-l) 

(4.13) 
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where we have defined Pn+i ^ P + Q, so M^+i = s, and Pi = p', so = 0. 

The dependence of the phase space on 1 — a; can now be traced by performing the change of 
variables 

= T7^'^ Mf = szn ■ ■ ■ Zi] i = 2,...n. (4.14) 
All Zi range between and 1, and 

dM^ . . . dMl = s'^-^z'^-'^zlzl ...zsdzn... dz2. (4.15) 
The two-body phase space Eq. (|A.9| ) is then 

dMP^+l■, Pi) = N{e) (zn z„-i . . . ^.+i)"^ (1 - z,y~^' dil,, (4.16) 

where dfli is the angular integral in the center-of-mass frame of the {ki, Pi) system, and A^(e) = 
1 

2(477)2-2-! • 

Hence, the n-body phase space Eq. ( [4.13|) can be rewritten as 



bn+iip + q;ki,..., kn,p') = 2n 



27r 

/' dZn 4"-'^-^"-'^^ (1 - Zn)'-'' ... f dZ2 Z^ (1 " Z2)'-'' ■ (4.17) 

Jo Jo 



N{e) 



„n— 1— ne. 



The dependence of the phase space on 1 — x comes entirely from the prefactor of 
the dependence on x and has been entirely removed from the integration range. Now, the 
amplitude whose square modulus is integrated with the phase space Eq. ( |4.17| ) is in general a 
function y4„+i(Q^, s; Z2, ■ ■ ■ , z^, ^1, ■ ■ ■ , ^n)'- it depends on 3{n + 3) — 10 + 1 = 3n independent 
variables for a process with two initial-state and n + 1 final-state particles, in which one of the 
masses (the virtuality Q"^) is not fixed. In the x ^ 1 limit, s — and the dominant contribution 
is given by terms which are most singular as s vanishes. Because of the cancellation of infrared 
singularities iTll, \An+ stronger singularity would lead to powerlike infrared 



■»0 



divergences. Less singular terms instead correspond to contributions whose Mellin transform 
vanishes in the large- limit. 

Hence only terms in the square amplitude which behave as s-"+o(<:) contribute in the x — >■ 1 
limit. In d dimensions, these terms pick up an s"""^ prefactor from the phase space Eq. (|4.17| ). In 
order to establish the result Eq. ( [4.4|) we must study the dependence of these terms on s*^. The 
case of a tree- level amplitude, i.e., purely real soft emission is in fact straightforward. Indeed, 
in this case, the amplitude manifestly does not contain any power of s"^, so the result Eq. ( [4.4| ) 
follows immediately, with a = 1. Notice that in this case, since the factors of s~'^ come from 
the phase space only, the sum over k in Eq. ([4.4|) reduces to the single term k = n: the power 
of [Q^(l — x)] coincides with the number of z integrations, i.e., with the number of emitted 
partons. 

To understand the physical meaning of this result, note that in the center-of-mass frame of 
the ki. Pi system the energy of the i-th of the n emitted partons is (see Appendix ^ 

. Mi: I ( Mf \ , , 

.. = .« = ^(l_^j. (4.18) 
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Now, concentrate on the emission of the n-th parton, whose energy is = ^{1 — Zn) = 
^ ( 1 — — ) . The singularity arises if the squared amphtude behaves as 1/(A;°)^. This behaviour 



is seen to be present e.g. for gluon emission in the eikonal approximation. In such case we get 



I A 



n+1 1 



M+1 



dk^ 



1 

2^ 



fco 



l+2e 



X] 



1 

(T 



Ax 



(4.19) 
(4.20) 



where the maximum value of /c° is attained when 



and is 



k^ 




(4.21) 



The X dependence of the result is thus simply determined kinematically in terms of the upper 
limit of the logarithmic integration. Eq. ( [4. 17] ) then shows that, thanks to the structure of 
the n-body phase space, this remains true in the general case of n-parton emission. 

We will now study how this result is modified by the inclusion of loops. To this purpose, we 
notice that a generic amplitude with loops can be viewed as a tree-level amplitude formed with 
proper vertices. Powers of s*^ can only arise from loop integrations in the proper vertices. We 
thus consider only purely scalar loop integrals, since numerators of fermion or vector propagators 
and vertex factors cannot induce any dependence on s^. 

Let us therefore consider an arbitrary proper diagram G in a massless scalar theory with E 



external lines, / internal lines and V vertices. It can be shown [|T^ that, denoting with P the set 
of E external momenta and Pe the set of independent invariants the corresponding amplitude 
Ag{Pe) has the form 



AciPE) = K{2nY5^''\Y.P)MPE) 

ii~L{d-i) j_ r ^1 



(4.22) 



AciPE) 



:r{I-dL/2)l[ f d(3i 
1=1 



[Pam^iL+m-i [DG(/?,Pi,)]'-'^/" 



In Eq. ( 4.23 ), Pi are the usual Feynman parameters, Pg{P) is a homogeneous polynomial of 
degree L in the Pi, Dg{P,Pe) is a homogeneous polynomial of degree L + 1 in the Pi with 
coefficients which are linear functions of the scalar products of the set Pe, i.e. with dimensions 
of (mass)^, and K collects all overall factors, such as couplings and symmetry factors. 

The amplitude Ag{Pe) Eq- (fl-23|) depends on s only through Dg{P, Pe), which, in turn, is 
linear in s. We can then determine in general the dependence of Ag{Pe) by considering two 
possible cases. The first possibility is that Dg{P, Pe) is independent of all invariants except s, 
i.e. Dg{P,Pe) = sdciP). In such case, Ag{Pe) depends on s as 



Ag{Pe 



I-2L+Le 



(4.23) 



where is a numerical constant, obtained performing the Feynman parameter integrals. The 
second possibility is that Dg{P,Pe) depends on some of the other invariants. In such case. 
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Ac{Pe) is manifestly an analytic function of s at s = 0, and thus it can be expanded in Taylor 
series around s = 0, with coefficients which depend on the other invariants. In the former case, 
Eq. (|4.23|) implies that the s dependence induced by loop integrations in the square amplitude is 
given by integer powers of s~^ In the latter case, the s dependence induced by loop integrations 
in the square amplitude is given by integer positive powers of s. 



Therefore, we conclude that Eq. (|4.4|) holds in general even when loop integrations are present. 
Unlike in the case of tree-level diagrams, however, the overall power of is no longer determined 
by the phase space only, and thus no longer equal to the number of emitted partons, from which 
it can differ by an integer amount. Specifically, each loop integration can carry at most a factor 
of Hence, at order the highest power of is n, given that either an extra real emission 
or an extra loop can give at most an extra factor but otherwise contributions proportional to 
all integer powers of s""^ up to n are allowed. Because in d dimensions the bare coupling constant 
tto has mass dimensions 2e, the dependence on Q^, which is the only other dimensionful variable, 
is fixed by dimensional analysis, and Eq. (|4.4| ) immediately follows. 

Furthermore, it is clear (compare Eq. ( [4.20| )) that each logarithmic Zi integration produces 
a 1/e pole. Since each angular integration gives a 1/e pole from the coUinear region, the highest 
soft singularity due to real emission in the coefficient C^^jj /T{—ake) in Eq. ( |4.4| ) is an e pole of 
order 2k — 1. Hence, C^^ has a soft singularity from the purely real emission term 

the sum Eq. ([4.3|). Because infrared singularities must cancel |Tl|] in the sum of all other 



terms C^°^ in the sum also has an order 2n pole in e. 

The whole argument can be reproduced for Drell-Yan-like processes with minor modifications 
which account for the different kinematics. Specifically, in the case of Drell-Yan we have massless 
partons in the initial state with momenta p, p' , while the Drell-Yan pair in the final state has 
momentum Q, so x = Q'^/s, with s = {p + p'Y- The recursive expression for the n-body phase 
space Eq. (|A.5|) now becomes 



rs 

d(j)n+i{p + p';Q,ki,...,kn)= / —^d4>2{p + p';kn,Pn)d(f)n{Pn;Q,ki,...,kn-i), (4.24) 
leading to 



d(j)n+l{Pn+i;Q,kl,...,kn) = / — C?02(P + p'; A;„, P„) / ""^ d(t)2{Pn\kn-l,Pn-l) 

Jq2 27r JQ2 271 

... / — ^d02(P3;A;2,P2)#2(P2;A;i,Pi), (4.25) 

where now Pn+i ^ P + p', so M^^^ = s, and Pi = Q. 

The change of variables which separates off the dependence on 1 — x is now 



^^ = T7i T^! Ml -Q' = {s-Q')zn...Zi; t = 2,...n (4.26) 



so all Zj still range between and 1 and 

dM^... dMl = {s- Q^y-^zl^^zlzl ...z^dzn... dz2. (4.27) 
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The two-body phase space Eq. 



IS now 

2 \-l+e 



{Ml, - Q') - (M^ - 



l-2e 



(4.28) 



where in the last step we have replaced {M^,) ^'^'^ by {Q^) ^'^'^ in the x —>■ 1 limit (compare 



Eq. (plD). 

We finally get 



{p + q;ki,...,k 

niP ) 



2n 



2tx 



iQ 



2N-n(l-e)/ 



Q 



2\2n-l-2ne 



dQn ■ ■ ■ d^li 



2e 



dZ2 Z2 



l-2e/ 



Z2 



sl-2t 



The dependence on 1 — a; is now entirely contained in the prefactor 

^1— 2ri+2ne 



(Q2)-"(l-)(s_g2)2n-l-2n.) 

Again, this proves the result Eq. 



i2ln— ne 



(4.29) 



(4.30) 



for real (tree-level) emission, since the corresponding 
amplitude is free of factors of (s — Q^)"*^. Note however that now a = 2. The physical interpre- 
tation is the same as in the deep-inelastic case, except that the upper limit of the logarithmic 
integration Eq. ( |4.19|) for the i-th parton is now attained when Mf = and is 



k' 



2 I s 



lQ'{l-x) 
Ax 



(4.31) 



Finally, loops can be included as in the case of deep-inelastic processes, with the possible cut at 
X = 1 now being given by a factor of (s — Q'^)^~'^^~^^'^. 

We conclude that Eq. (^4.4| ) holds for both deep-inelastic and Drell-Yan-like processes, with 
a = 1 in the former case and a = 2 in the latter. This difference is merely a reflection of the 
underlying kinematics, Eqs. ( [4.21|) and (|4.31| ), but in both cases the soft logs {i.e. powers of 
(1 — x)""^) can be traced to the fact that the phase space variables Mf can all be related to a 
single dimensionful variable s (DIS) or s — (DY) through dimensionless variables Zi which 
range from to 1, according to Eqs. ( [4.14| , 06D . 

It is interesting to ask which are the kinematic configurations which contribute to the x 1 
limit: in fact, they turn out not no be the same in deep-inelastic and Drell-Yan. Indeed, consider 
a contribution to the cross section that involves n massless partons with momenta fci, . . . , /c„ in 
the final state. We have 

p + p' = q + ki + ... + kn. (4.32) 
Squaring both sides of Eq. ( [4. 32] ) we get 



s(l - x) = ^ fci ■ fcj + 2 ^ g ■ ki. 

i,j=l i=l 



(4.33) 
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The quantities 

q-ki = (go - ^qi " Q' cos 0,^ (4.34) 

are positive semi-definite, and vanish only for = 0. Hence, for Drell-Yan processes, when x 
approaches 1, all fc^'s must go to zero. i.e. all emitted partons must be soft. 
On the other hand, in the case of DIS we have 

p + q = ki + ... + kn (4.35) 

(where ki = p'), which gives 

^ ^ ' = T.ki-k,= Y: k^k^ (1 - cos %), (4.36) 

^ i,j=i i,j=i 

where 6ij is the angle between ki and kj. Now take the limit x — 1. The most general solution 
to Eq. (^4.36|) is not ki=0 for all i, but rather 

k^ = for l<^<n 

6ij = 0;k^,k'j arbitrary for n < i, j < n. \ ■ J 

In other words, for DIS the kinematical configurations that correspond to a; = 1 include not only 
the case when all partons in the final state have small energies, but also the cases when a subset 
of them are collinear to each other. 



4.2 Renormalization group improvement 

We now study the restrictions that renormalization-group invariance imposes on the cross sec- 
tion. As already mentioned, our only assumption is the standard factorization Eq. ( |2.2| ), namely, 
that C{N,Q'^/fi'^,as{fi'^)) can be multiplicatively renormalized. This means that all divergences 
can be removed from the bare coefficient C^^^ {N, Q'^,ao, e) Eq. ( [4.2|) by defining a renormalized 
coupling as{fi'^) according to the implicit equation 

«o(/i', «.(/i'), e) = /x'^ a.(/x') e) (4-38) 

and a renormalized coefficient 

C (iV, Q^fi^ «,(/i2), e) = Z(^)(iV, «,(/i2), e) C(°) (iV, Q\ «o, e) , (4.39) 

where fi is the renormalization scale and ^''^"^^(^^(/i^), e) and Z^'^^N, asifJ^"^), e) are computable 
in perturbation theory and have multiple poles at e = 0. The renormalized coefficient func- 
tion C (A^, Q^Z/i^, as(/i2), e) is finite at e = 0, and, because the renormalized coupling as is 
dimensionless, it can only depend on through Q^Z/i^. 

In a mass-independent renormalization (and factorization) scheme such as MS, both Z^""^ 
and Z^'-^'^ are series in as with Q^-independent coefficients. Therefore, the physical anomalous 
dimension is 

7 = In C {n, gV/i', e) = -e«o/- In C^'^ {n, Q\ Cq, e) , (4.40) 
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where we have exploited the fact that, for dimensional reasons, C^^^ can only depend on 
through the combination aoQ"^*^. Q Equation ( |4.40| ) implies that in d dimensions the physical 
anomalous dimension 7 viewed as a function of ao admits an expansion in powers of aoQ~'^'^, 
while as a function of a<i(/i^) it admits an expansion in powers of (Q^Z/i^)"*^- 

In the large limit, the dependence of the bare coefficient function on is given by Eq. ( [4.3|) . 
But Eq. ( 4.40 ) implies that 7 has the same property, namely 



7 



ao,e 



i=l j=0 



EE 7., (6) [{QT'c^o 

i=l j=0 



+ 



ao 



(4,41) 



Hence, in the large N limit the d-dimensional physical anomalous dimension is a power series in 
Q~'^'^ao and {Q'^ /N'^)~''ao, with A^-independent coefficients. 

The renormalized expression of the physical anomalous dimension is found expressing in 
Eq. ( [4.41| ) the bare coupling in terms of the renormalized one by means of Eq. ( |4.3(j| ): this yields 
an expression of the physical anomalous dimension which is finite in the limit e — > 0. Now, the 
function 



is manifestly renormalization-group invariant: 

2daQ 



(4.42) 



/i 



0. 



It follows that 

«o(QV/i',«.(/^'),e) = «o(l,«s(Q'),e) = «,(g2)z("»)(a,(g2)^. 
Explicitly, noting that the d dimensional beta function is given by 
o d 



(4.43) 



(4.44) 



= -ea,(^2) -/5oa^(/i2) +0(af) 



(4.45) 



^Note that the standard anomalous dimension 7^^ Eq. ( |2.4| ) is instead given by 



7^P(^,a,(A*^)) = -: 



d In fi^ 



due to the fact that the bare coefRcient C*^°^ is independent of /i^. It follows that 7^^ = 7 if and only if 
Z^°'''\as{fM^)T e) — 1, i.e. if the four-dimensional beta function vanishes. 
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in terms of the four- dimensional beta function of Eq. ( p.l3| ), and exploiting the //-independence 
of ao, Eq. ( |4.38| ) gives 



exp 



" da' P{a') 



a' ea' — (3{a') 



1 + 



/3oas 



+ 



(4.46) 



The renormalized anomalous dimension is found by rewriting Eq. ( [4.41| ) as 



i=l j=0 



ao(l,as(QVA^"),e) 



ao(l,tts(Q^),e) 

and then re-expanding in powers of the renormalized coupling: 

oo m 

7(iV,QV/i',«s(/x2),e) = E E7L(e)«r"(Q')<(QViV'^)- 



(4.47) 



(4.48) 



m=l n=0 



The physical anomalous dimension is thus seen to be a power series in ^^(Q^) and aglQ"^ /N"-), 
with coefficients determined by 7ij(e) and the renormalization constant Z^°'''\ If as{Q'^/N'^) 
is re-expanded in terms of a<((Q^), the physical anomalous dimension is seen to be of order 
tts(Q^) (and thus vanish when as(Q^) — > 0), as it ought to. 

However, we cannot yet conclude that the four- dimensional physical anomalous dimension 
admits an expansion of the form of Eq. ( |4.4(j| ), because the coefficients 7^n(^) are not necessarily 
finite as e ^ 0. In order to understand this, it is convenient to separate off the A^-independent 
terms in the renormalized anomalous dimension, i.e. the terms with = in the internal sum 
in Eq. ( |4.48| ). Namely, we write 

7 (iV, QV/i^ a.(/x2), e) = 7^'^ («s(Q'), as{Q'/N^), e) + 7^^^ («.(g'), e) , (4.49) 
where we have defined 



f') (a,(Q2),a,(QVA^"), 



EE7.,(e) [ao(l,as(Q'),e)]'"' [ao(l, a.(QVA^"), e) 
i=i j=i 

00 00 

m=0 n=l 



and 



r- («.(Q^),e) 7.0(e) 

i=l 



ao(l,«.(Q2),e)l^= E7l(e)«r(Q'). 



(4.50) 



(4.51) 



m=l 



Whereas 7 (iV, a^, e) is finite in the limit e — > 0, where it coincides with the four- 

dimensional physical anomalous dimension of Eq. (|2.3|) , 7^'^ and 7*-'^^ are not necessarily separately 
finite as e — i> 0. However, Eq. ( [4. 491 ) implies that 7''''' and 7*^^^^ can be made finite by adding and 
subtracting to them a suitable counterterm Z^^^ (a;s(Q^),e): 



7^') {as{Q'),as{QVN^), 
1^^ («s(Q'), 



7^') (a.(g'), a,(gViV'^), e) - Z^^) (a^m, e) , (4.52) 
7(^)(a,(Q2),6)+zW(a,(Q2),6) (4.53) 
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where the functions 7*^'^ and 7*^^^^ have a finite e — hmit. The counterterm cannot depend on 
N because of the A^-independence of 7'-'^^. A possible choice is therefore 

{a^m, e) = 7^'^ e) . (4.54) 

With this choice, 7''''' is clearly finite for N = 1, where it vanishes; because of the A^-independence 
of the counterterm, 7'-') will then be finite for any A^. Other forms of the counterterm can differ 
from this only in the choice of the finite part, i.e., can be obtained by adding to Z^'^^ Eq. ([4.54|) 
a finite function of as{Q'^). Therefore, they correspond to the possibility of reshuffling finite 
A^-independent terms between 7*^'-* and 7*-^-'. The choice Eq. ( [4.54| ) is characterized by the fact 
that with this choice 7« is 'purely logarithmic', in that it vanishes at A^ = 1 where InA^ = 0. 
We thus get 

7^') {a^m, c^siQ'/N'^), e) = 7^'^ («s(Q'), a.iQ'/N^), e) - 7^') (^.(Q^), ^.(Q^), 6)(4.55) 
7^^) {a,m, e) = 7^^^ («s(Q'), e) + 7^'^ («s(Q'), «s(Q'), e) , (4.56) 
and the physical anomalous dimension is 

7 (at, QV/i', «.(/i'), e) = 7^'^ («s(Q'), asiQVN'^), e) + 7(^H«.(Q'), e) + 0{1/N) 

= 7^'^ («.(Q'), a.iQVN'^), e) + OiN^), (4.57) 

where now both 7^'^ and 7^^^-' are finite as e — 0, and 7''''' provides an expression of the resummed 
physical anomalous dimension in the large A^ limit, up to non- logarithmic terms. 

It is apparent from Eq. ( |4.56| ) and the definition Eq. ( |4.51| ) that 'y^^^ is a power series in 
tts(Q^) with finite coefficients. In order to understand the perturbative structure of 7*^'-' as well, 
we notice that 

7(')(«.(Q'),a.(QVAr'^),e) = /"^^ ^ ^? (^.(g^), a,(gVn), e) (4.58) 



where 



d 



/3^'n«.(/i'),e)^^^7^') («s(g'),«.(/i^),e) , (4.59) 



d 



where /3^°'^(as) is the dimensional beta function of Eq. ( ^4.45|) . It immediately follows from 
Eqs. (|4. 50144.591) that g is a power series in ^^(Q^) and ^^(/i^): 

00 00 

g («,(Q2), «,(/i2), e) = E E w(e) aT{Q') <(/x'). (4.60) 

m=0 n=l 

The sum over powers of as(yU^) in Eq. ([4.60|) starts at n = 1 because the expansion of 7*^'^ 
contains at least one power of as{Q'^/N). The perturbative structure of 7*^'^ is thus dictated by 
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Eqs. ( 4.58 , 4.60 ). Because of Eq. ( 4.57 ), this provides us with the desired all-order generalization 
of the next-to- leading log resummation formula Eq. ( |2.23| ) discussed in Sect. ^ 

Our final result for the four- dimensional resummed physical anomalous dimension is thus 

rN°- fin 

7'-^^(iV, as{k^)) = / —g{as{e), as{e/n)) + 0(iV°) (4.61) 
Ji n 

oo oo 

g{as{e). asik'/n)) = E E W«r(fc')<(fcV^); 9n = 0. (4.62) 

771=0 71=1 

Exploiting the all-order relation between leading In jj: and leading ln(l — x) resummation derived 
in Sect. |^, Eqs. ( p.l9| , p.20D , this result can be equivalently cast as the Mellin transform of a 
resummed splitting function P^^^(a;, as(fc^)): 



Y^'iN^asik')) = / dxx 
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^-ip'^'^^(x,«,(A;2)) + 0(iV°) 
g{as{e),a,{k'il-xr)y 



1 — X 



f.^^^^-^9io.s{k%a,{l^')). 



p=0 



(4.63) 
(4.64) 

(4.65) 



These resummations generalize to all orders the next-to-leading log results Eq. ( |2.23| ) and 
Eq. (|2.11| ), respectively. The two forms Eq. (|4.61| ) and (|4.63|) of the resummation are equivalent 
up to non-logarithmic terms, and can be used to compute the resummed evolution factor Eq. ( |2.6D 

K{N; Ql Q') = exp J^^ -^^'iN, ^.(fc^)). (4.66) 

Resummed expressions for the coefficient function and anomalous dimension in any factorization 
scheme can be obtained from Eq. (f4.66| ), as we will discuss in Sect |^. 



4.3 Next-to-leading resummation 

Let us now study our all-order generalized resummation at the leading and next-to-leading level. 
This will also help us to understand the relation of our resummed result to the stronger one 
proven in Ref. The lowest-order term in Eq. ( |4.60| ), used in the expression of the resummed 
anomalous dimension Eqs. ( [4.571 - |4l5^ ), leads to the leading log resummation: 
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(0 



-fl'oi 



da 
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goi , (3oasiQy N^) + e /n^/Ar^^ 



(4.67) 



which reduces to the leading-log result Eq. ( |2.28| ) as e ^ 0. 
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Recalling that 7^'' Eq. (|4.55|) is independent of as(Q^) and is free of constant {i.e. a^- 
independent) terms, Equation ( [4.67| ) uniquely determines the leading log contribution to 7^'^: 



/^^/WOViV^V (4.68) 

This shows explicitly that 7''''' is an analytic function of as{Q'^ /N), but does not have a finite 
e ^ limit. By contrast, 71''' Eq. (|4.671 ) does have a finite e — >■ limit, but the limit is no longer 
an analytic function of both as{Q'^/N) and asiQ"^). 

Subleading resummed terms are found by integrating according to Eq. (|4.58|) higher order 
contributions in the expansion Eq. ( |4.6CI| ), and consistently including higher-order contributions 
to the beta function. In general, the term gij will generate a next*"'"-' ~ ^-to-leading order resumma- 
tion (when the beta function is included at leading order). The NLO resummation of Eq. ( |2.23| ) 
is recovered if the only subleading contribution is the go2 term in the expansion Eq. ( [4.60| ), i.e., 
if ^711 = 0. 

As discussed in the introduction, an all-order generalization of the resummation Eq. (|2.23|) 
was derived in Ref. p. This resummation has the form of Eqs. ( |4.57| - ^^58|) , but with g a function 



of Q;s(/i^) only, i.e., with all gmn = when m > in Eq. ( [4.60| ). In such case, knowledge of 

the first k coefficients gon, 1 < n < k fully determines the next^~^-to- leading log resummation. 

These coefficients, in turn, can be determined by performing a next^~ ^-to-leading fixed-order 

computation. The more general resummation of Eq. (|4.60| ) is less predictive. Indeed, if all 

coefficients gmn are nonzero, knowledge of the k{k+l)/2 coefficients with m+n < kis necessary in 

h—i fc(fc+i) , 

order to determine the next -to-leading log resummation. In turn, this requires a next 2 

to-leading fixed-order computation. So, while the leading-log resummation is determined by the 

0{as) result, an 0{al) computation is required in order to determine both coefficients (720 and 

gii which control the next-to-leading log resummation, and so forth. 

The conditions under which the more restrictive result of Ref. P] holds can be understood 

by comparing to our approach the derivation of that result. The approach of Ref. is based 

on assuming the validity of a factorization formula which is more restrictive than the standard 

collinear factorization Eq. ( |2.2| ). This factorization was proven for a wide class of processes in 

Ref. and it implies that the perturbative coefficient Eq. (|2.2| ) in the large- limit can be 

factored as: 

C{N,Q'/fi\a,{fi')) = C«(gV(/i2A^'^),a,(/i2))C(^)(QV/i^a,(/x2)). (4.69) 
or equivalently 

C^'\N,Q',ao) = C^'^'\Q^/N'^,ao,e)C^''''\Q\ao,e). (4.70) 
Then, the physical anomalous dimension Eq. ( ^.401 ) manifestly becomes 

7 {n, gV/i', «s(/i'), e) = 7^') (gV(iVV), e) + 7^^^ {QV «s(/x^), e) , (4.71) 

i.e., 'j^^'' depends on asiQ"^ /N"') only. Proceeding as above, one then ends up with the resum- 
mation Eqs. (|4.57| - [4.58|) , but with g = g{as{ii'^), e), or, equivalently, all gmn = when m > 0, as 
advertised. 
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By running this argument in reverse, the factorization Eq. ( 4.69 ) is thus seen to be a necessary 
and sufficient condition for the vahdity of the more restrictive factorization formula of Ref . p[ , 
where g = g{as{n'^))- This factorization is in general rather nontrivial. Indeed, recalling the 
expansion Eqs. ( ^I2| , [4.3| ), Eq. ( ^4.7CI| ) is seen to be satisfied if and only if in the large- iV limit the 
coefficients C^\e) can be written as 



C^'\^) = j:Fk{e)G„.^,{e) (Q^ 



k=0 



-ke 



(4.72) 



Whereas an all-order proof of Eq. ( 4.72|) is highly nontrivial, the structure of our more general, 
but weaker resummation Eqs. ( f4.57| - |4758| ) implies that a fixed order proof of the factorization is 
in fact sufficient for the corresponding resummed result to hold. 

Consider the simplest nontrivial case, namely, the next-to-leading log resummation: this is 
determined by knowledge of the coefficients (720 and Qu. The stronger resummation holds when 
5^11 = 0. But in order to check whether this is the case or not it is sufficient to verify that the 
desired factorization Eq. (|4.72| ) holds to lowest nontrivial order, namely that 



Ci'\N,Q\e) = G,ie){Q 



2\-2e 



Q 



~2e 



Fi(e)Gi(e)(g^)- ^ + F,{e) ^ + 0{1/N). (4.73) 
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Equivalently, we can compute the coefficients goi, g2o and gn which determine the next- 
to-leading log resummation by comparison with fixed-order results: in order to calculate three 
coefficients, we need a next-to-next-to leading fixed order computation. Indeed, expanding 
Eqs. (|4.61| ) in powers of as{Q'^) to order we get 



= -fl'oi«s(Q^) + 



(4.74) 



^ In^ - - i9n 
2 ^ 



■ fei J 3 1 
3 



I + Po9o2 + 



Om) In — 



Pi9o2 Plgm ^ 



Pi 



In- 



A^ 



+0{at) + 0{N'). 



Hence, after having determined g^i from the leading log result, the O(a^), O(lnA^) term only 
determines the combination gn + (702, but the O(a^), 0(ln^ A^) provides an independent hnear 
combination. 

The physical anomalous dimension is in turn determined by the coefficient function C 
Eq. ( |2.2| ) and the standard anomalous dimension 7"^^ Eq. ( |2.4j) according to Eq. ( |2.5|) , which 
through order gives 

(4.75) 

where 7^^^ and Ci are the order-a* coefficients in the expansion of 7^^ and C respectively. 

The coefficients C are known up to 0{al) |]l5l for Drell-Yan |]l3l and deep-inelastic [l^ 
processes, while 'j^^ is known (in the MS scheme) up to O^a^) [|l^, and has been very recently 
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determined to O(a^) in the large limit |17|. However, knowledge of the O(a^) anomalous 
dimension is in fact sufficient to determine gu, thanks to the fact that in the MS scheme the 
anomalous dimension is linear in IniV to all orders |T^, so its knowledge to ©(a^) is sufficient 
to establish the vanishing of gu. 

To see how this works, consider explicitly the case of the deep-inelastic structure function 
F2. The coefficients in the MS scheme are 



CI 




3/o-4C(2)-9]+0(l/iV) 

4/3 - 9/2 
11 /367 



27' 
16C(2) + 



4C(2) h 



^ + 18C(2)-4C(3))/o' 

3155 22 , , 
-^ + -C(2) + 20C(3) 

|+0(iV°), 



(4.76) 



(4.77) 



where ( is the Riemann ( function, Cp and Ca are the usual quadratic Casimir operators, and Ip 
have been defined in Eq. ( |3.2| ) and, as proven in Sect. ^ contain the logarithmic N dependence. 
Note that 0{N^) terms have been included in Eq. ( [4.76|) in view of the fact that Cf appears 
in Eq. ( |4.75| ). The gluon coefficient to next-to-leading order is 0{1/N). The MS anomalous 
dimension in turn is 



goAP 
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/o + 0(ArO 
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UfCp + 
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2C(2) CaCp 



/o + 0(ArO 



(4.78) 
(4.79) 



while 7''^ and 7^^ are 0{1/N). The gluon anomalous dimension 7^^ does contain soft logs, but its 
contribution to F2 only goes through via interference with C^, and it is thus power suppressed. 
Therefore, we see that through next-to-leading order there is no operator mixing, and only the 



quark sector contributes to the large limit. Substituting in Eq. ( [4. 751) , and solving for gu and 
(7o2 we get, after tedious but straightforward algebra 
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= 0. 



(247s + 58)Cirn/ + {12n^ - U2-fE - 367)CaCf 



(4.80) 

(4.81) 
(4.82) 



In Drell-Yan there is likewise no mixing, and the check that gu = 0, performed in Ref. 
follows the same lines. 

Hence, by combining available fixed-order results with our general resummation Eq. ( [4.61| ), 
we reproduce the next-to-leading log resummation Eq. (|2.23| ). Beyond next-to- leading order 
the general result Eq. (|4.61|) holds, and will have to be used for the full matrix of two-by-two 
resummed physical anomalous dimensions. 
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5 Momentum-space resummation 



We have seen that resummed results can be compactly expressed in terms of the physical anoma- 
lous dimension Eq. ( |4.61| ), or in terms of a resummed splitting function Eq. ( |4.64| ) obtained from 
it by inverse Mellin transformation, up to non-logarithmic terms. In the former case, the resum- 
mation is cast in Mellin space, and one resums large In (at large — > oo). In the latter case, 
the splitting function is defined in momentum space, i.e., as a function of the momentum fraction 
X and one resums large ln(l— a;) (at large x — > 1). In this section we address the question whether 
it is possible to construct a resummation of ln(l — x) directly at the level of the momentum-space 
cross section a{x, Q"^). Indeed, even when the resummation is performed at the level of splitting 
functions, the cross section is written as the inverse Mellin transform of the exponentiated result 
Eq. ( ^.111 ). It would seem more natural if the momentum-space cross section (t(x, Q^) could be 
written directly as the exponential of a suitable leading, next-to-leading,. . . ln(l — x) function. 

A general way of performing the inverse Mellin transform of a leading-, next-to-leading, . . . , 
In function up to terms which are suppressed by powers of 1— x was given in Sect. |], Eq. ( p.25| ). 
Using this result, it is easy to rewrite the A^-space resummation factor Eq. ([4.66|) as the Mellin 
transform of an x space resummed expression: 



tiki Jo in-l)\ 



1 



1 -xrfln'^d -X 




t^7«(«.(^^),«s(fc^(l-a:r)/ 
Ql 



(5.1) 



+ 0(ArO), 



where 7« is the resummed physical anomalous dimension Eq. ( [4.58| ), and, even though the 
resummation factor K^^^ is defined up to non-logarithmic terms, we have normalized it by the 
requirement that K^°^ ^ 1 as — 0. This immediately leads to the identification of the 
momentum-space resummation factor: 
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(5.2) 



We can now derive a closed form of the sum of the exponential expansion in Eq. (|5.2|) , using 
a representation of the -|- distribution given in Ref. 0]: 
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(5.3) 
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where 6 is the Heaviside function. The equahty holds in the sense of distributions: it is proven 
by folding K^'^^ with a test function, integrating by parts, and noting that 0{—ri) = and 
7*^'-*(as(Q^), = because of Eq. ([4.55|) . It is also convenient to define a momentum- 

space resummation exponent E^'^^{x;Ql,Q'^) through the implicit equation 



lim — ^ 

r7-»0+ ax 



9{l-x-r])expE'''%x;QlQ' 



(5.4) 



Terms with n = 1, 2, ... in Eqs. (|5.2| , ^.3|) generate leading, next-to-leading,. . . ln(l — x) contri- 
butions to the resummation factor. Hence, if 7^^^ is computed at the leading log level Eq. ( p^.28D , 
only the n = 1 terms should be retained. 
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(5.5) 
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(5.6) 

At next-to-leading and higher orders, the action of subsequent terms in the series of derivatives 
Eq. ( |5.3| ) should be accordingly truncated so that the resummation exponent E^'^^ is computed 
to the given order. 

We would thus be led to conclude that Eq. ( [5.5| ), and its obvious subleading generalizations, 
provide us with the desired form of the leading-, next-to- leading-,. . . , ln(l — x) resummation 
factor, which can then be used to construct the corresponding resummed momentum-space cross 
sections. Surprisingly, however, it turns out that the formal expression Eq. ( ^.5|) is actually ill- 
defined: the corresponding resummed cross section diverges. This difficulty was first pointed 
out in the fixed coupling limit in Ref. [0. Also surprisingly, this divergence is a consequence of 
the leading- (next-to- leading,. .. ) ln(l — x) truncation: it is absent if the Mellin transform of 
the leading- (next-to- leading,. .. ) In-^: result is performed up to power suppressed terms, i.e., 
by retaining all terms in the sum over n Eq. ( p.3|) . 

To see this divergence, recall that the resummed evolution factor K^^^ according to Eq. ( |2.6| ) 
satisfies 



a{z,Q')= f'^K^^^{x-Ql,Q')a(-,Qi) ^ f' dx K^^^x; Ql, Q') t{x) 

Jz X \X J Jo 



where 
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r{x) = -a{-,Qt\ e{x-z). 
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Integrating by parts and noting that r(0) = 0, we get 







'rfa;^^ expE'^^^(x;g2,g') 



dx 



(5.7) 



(5.8) 



(5.9) 



26 



Specializing to the leading log resummation Eq. ( ^.5|) 
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where we have expanded out the leading log exponent Eq. ( |5.6| ) in powers of a^, and 0(a;^) 
denotes the next-to-next-to leading term in this expansion, which is in fact of order {a1{Q'^) — 



X] 



ali^Q"^) ln^(l - x)"" [1 + 0(a,(Q2))]. Explicitly, the coefficients of the first 



two orders of this expansion are Go = fi'i/Z^o and Gi = —gi/2. 

Now, as X ^ 1, the 0{as) contribution to E^^ diverges as ln^(l — x), so exp E'^'^ diverges 
faster than any power of 1 — x. However, r(x) is manifestly regular as x — 1, so the integral 
Eq. (|5.9|) does not exist. In fact, matching with the leading-order anomalous dimension shows 
that gi = —Cp/'rc < 0, so the integral diverges for < Ql; however (regardless of the sign of 
5^1) the presence of this divergence means that the perturbative expansion of a in powers of 
has vanishing radius of convergence. 

This can also be seen explicitly by using the leading log K^^ Eq. ( |5.5|) in Eq. ( ^.7|) : 
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The generic term in the sum over p is 
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a,(Q2) \2p+l 



where in we have Taylor expanded r(x) about x = 1, and neglected terms which are suppressed 
by powers of 1 — x, consistent with the fact that the resummed distribution K^*^^ is determined 
in the same approximation. The factorial growth of Cp Eq. (|5.12 ) implies that the sum over p 
in Eq. ( |5.11|) diverges factorially. 

Higher order contributions in the expansion of K^^ in powers of «s will make this divergence 
worse. Indeed, the 0(aj) contribution is proportional to ln'^(l — x), so the factorial divergence 
of the corresponding expansion is more severe. Furthermore, the expansion of K^^ amounts 
essentially to the expansion of ln(l + a<j/3oln(l — x)") in powers of as, which, since x < 1 is a 
series of terms with the same sign, so these divergences add up. 

In fact, the leading log exponent Eq. ( |5.(j| ) blows up when ln(l — x)*^ = ^^7^- This diver- 
gence corresponds to the Landau pole in the coupling constant as{k'^{l — x)") which appears 
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in the x-space resummation factor Eq. (|5.2| ). Because of the Landau pole, the expansion of 
the resummed exponent E in powers of has a finite radius of convergence. However, as first 
clarified in Ref. [^, the factorial divergence Eq. is entirely unrelated to the Landau pole. 

Indeed, we have seen that it is present even if we retain only one term in the expansion of the 
resummed exponent: rather, this divergence plagues the series expansion of the exponential, 
when integrated term by term. 

Now, we exploit the power of our all-order relation between leading ln(l — x) and leading In jj^ 
series to show that this factorial divergence is just a byproduct of the truncation to leading (or 
next-to-leading, etc.) level of the full expression Eq. ( ^.2| ) of the x-space resummation factor. 

To this purpose, let us again concentrate on the leading log expression, and separate off the 
first few terms in the expansion of the resummed exponent in powers of a^: retaining all terms 
in the sum over n Eq. (|5l^ ), we get 
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where in the last step we have used the definition Eq. (|3.13| ) of A*^"-*, and the limit ?7 ^ should 
be taken after performing all derivatives. 

Now, it is easy to compute the inner "leading order" sum (index p): exp acts as a finite 
translation operator. We get 
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We see therefore that, due to the inclusion of subleading corrections of the sum Eq. 
generic term in the sum over q is now 
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(5.16) 

77=0 



where we have again Taylor expanded r(x) about x = 1. 

Comparing this with the previous result Eq. ( ^.12| ) we see that the power in the denominator 
is replaced by a Gamma function. This is enough to tame the factorial divergence as we now 
show. We observe that 
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hence, Cq are the coefficients of the Taylor series expansion of an entire function, which, therefore, 
has infinite radius of convergence (note that if we had used the result Eq. ( [5. 12] ) we would have 
found instead a finite radius of convergence). This in turn implies that 



and consequently 



li^ 7%^7r = hm ^ %ti = 0, (5.18) 

lini ^^,4- = lim ^ %±i = 0. (5.19) 

This means that the coefficients C2q no longer diverge factorially: in fact, they are also coefficients 
of the Taylor series expansion of an entire function, with an infinite convergence radius. 

Proceeding in the same way, we can then show that the sum over r in Eq. (|5.14|) , i.e., the 
expansion of the expO(ag) contribution, is also convergent, and so forth. Therefore, we conclude 
that order by order in perturbation theory, the resummed result Eq. ( ^.2|) leads to a convergent 
perturbative expansion, but if this is turned into a leading-, next-to- leading,. .. In (1 — x), a 
meaningless divergent expansion is found. 

In summary, we have seen that the leading-, next-to-leading,. . . ln(l — x) expressions of 
the momentum-space cross section a{x, Q^) are all ill-defined. However, if the resummation is 
formulated at the leading-, next-to-leading,. . . In -i: level, and the resummed a{x, Q^) is defined as 
the inverse Mellin of the ensuing result, up to power corrections, then finite results are obtained. 



6 Factorization schemes and matching 

The resummed results derived in Sect. ^ provide us with a prediction for the resummed scale 
dependence of the physical cross section, expressed by the physical anomalous dimension 7. In 
practice, we would like to have resummed predictions for the coefficient function C Eq. ( p.2| ) and 
standard anomalous dimension 7"^^ Eq. (|2.4| ). Because of Eq. (|2.5| ), that relates 7 to 7"^^ and C, 
all the relevant information is contained in the physical anomalous dimension, and the resummed 
anomalous dimension and coefficient function can be determined from it. In order to do this, we 
must suitably combine resummed and unresummed results, and choose a factorization scheme. 
The simplest choice consists of choosing a physical factorization scheme, defined by 

aiN,Q') = FiN,Q') (6.1) 
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so that 



C(iV,l,a,(Q2)) 



7^^(iV,a,(Q^))=7(iV,«.(g')). 



(6.2) 



Since now the physical and standard anomalous dimension coincide, Eq. ( [4.6 1| ) gives us a re- 
summed expression of this anomalous dimension, up to non-logarithmic terms. A full expression 
of the anomalous dimension, valid both at large and small N, can be obtained by combining 
these resummed results with standard fixed-order ones, and subtracting the double counting. 
Specifically, a next-to-leading order and next-to-leading log expression of the physical anoma- 
lous dimension is given by 



7(iV,a,(A;2)) = a,{eW\N) + al{e)j^'\N) + 



^" dn 



1 n 



goia,{P/n) + go2al{k'^/n 



-goiasik') In ^ + aHP) l^ln" - go2 In ^ 
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where 7*^°-' is the standard leading-order anomalous dimension, 7*^^) is the two-loop physical 
anomalous dimension, obtained using Eq. ( |2.5| ) from the next-to-leading order anomalous dimen- 
sion and one-loop coefficient function, and the term in square brackets is the double-counting 
subtraction of the logarithmic contributions to 7^°-' and 7'-^^ 

A different resummed expression can be obtained by using the expression Eq. ( [4.63| ) of the 
resummed physical anomalous dimension in terms of the resummed splitting function Eq. ( [4.611 ). 
In this case, Eq. ([0|) is replaced by 
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where P*^*)(x) are defined through 
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and the relation between g and g was given in Eq. ( [4.65| ) (or in Eq. ( p.27[ ) to next-to-leading 
order). The two forms of the resummed anomalous dimension Eq. ([6.3[) and Eq. ([6.4|) coincide at 
order and at the next-to-leading log level, but differ by next-to-next-to leading terms, as well 
as by 0{1/N) terms which are at least 0(af). Hence, their difference might be used to estimate 
the impact of the next-to-next-to-leading resummation in a factorization-scheme independent 
way. 

Resummed coefficient function and anomalous dimension in a generic scheme can be obtained 
by generalizing the conventional factorization scheme change at the resummed level. A scheme 
change is performed by redefining the parton distribution by a finite function Z{N,as), which 
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can be viewed as a redefinition of tlie finite part of tlie multiplicative renormalization Z^'^'> of 
the hard coefficient, Eq. ( f4.39| ) 

F'{N,fi') = Z{N,as{fi'))F{N,fi^). (6.6) 
Upon the scheme change, the anomalous dimension and coefficient function change by 

In C'{N, Q/^i\ as{^^)) = In C(iV, Q/^i\ ^.(/i^)) _ In Z{N, as{^?)) (6.7) 

(6.8) 



7'^^(iV, as{^i'')) = 7^^(iV, asifJ.')) + /x^^ In Z{N, a,(/x^)). 
A conventional (unresummed) scheme change is performed by computing Z to fixed order in 



Z-(iV, «,(/i2)) = 1 + Zi(iV)«,(/i2) + Z,{N)al{fi^) + .... (6.9) 

The term Z^ affects the next'^-to-leading and higher order contributions to the anomalous dimen- 
sion. At the resummed next*^-to-leading log level, the physical anomalous dimension is computed 
including terms of 0(a""'"'^ In"" A^) to all orders in a.^- Such terms are reshuffled between the co- 
efficient function and the parton distribution by performing a scheme change Eq. ( |6.6|) with 



Z'-^^(Ar,«,)=exp 



zq as In 



AT" 



zi as In 



+ asZ2 as In 



Zk 



i=2 



C(s In -— 



(6.10) 
(6.11) 



where Zk affects the next'^-to-leading log and higher contributions to the resummed anomalous 
dimension and coefficient function C. 

Even though any resummed scheme change of the form Eq. (|6.10|) is in principle acceptable 
(i. e., it can be obtained by redefining the finite part of Z^"^^ Eq. ( [4.39D ), in practice it is convenient 
to consider in particular resummed scheme changes which preserve the structure Eq. ( |4.61|) of 
the resummed anomalous dimension, i.e. such that after scheme change Eq. ( |6.8|) the anomalous 
dimension can still be written in the form of Eq. ( ^.61| ), and the scheme change only affects the 
numerical values of the coefficients gij. Such scheme changes are performed by choosing Z to be 
of the form 

Z^'^{N,asifi'))=exp / — / ^ (z(a,(P), a,(A;Vn)) - z(«.(P), «.(P))) (6.12) 

00 00 

z{as{[i^), as{[i^ln)) = % "sl/^V"-); ^11 = 0. (6-13) 

i=0 j=l 

where A is an arbitrary fixed scale, and the subtraction is necessary in order to ensure that 
XwZ"^^^ be free of linear InA^ terms, in agreement with Eq. (|6.11|) . 

Resummed anomalous dimension and coefficient function in any scheme can in general be 
obtained from the full physical scheme anomalous dimension, which combines resummed and 
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unresummed results (such as Eg. ( |6.3|) ) by performing the combination of a resummed ( |6.10 ) and 
unresummed (|6.9|) scheme change, 



Z{N,a,) = Z\N,as)Z'''\N,a, 



(6.14) 



In practice, in the sequel we will take Z'^'^^ to be of the form Eq. ( |6.13| ). Then, Z^'^^ contains 
higher order powers of In but no linear terms in InA^ (nor, of course, non-logarithmic terms). 
Furthermore, we will take Z^ to contain at most a linear term in IniV, but not higher order 
powers. In such case, there is no overlap between a resummed next'^-to-leading log scheme 
change Z'^'^ and an unresummed next^-to-leading order scheme change Z": the unresummed 
and resummed scheme changes affect different contributions in an expansion of the anomalous 
dimension in powers of as and In A^. 

Let us now consider explicitly the commonly used MS scheme. This scheme is characterized 
by the fact that Z^'^'^ Eq. ( [4.39|) is a series of pure poles inl/e = l/e — 7^; + log(47r). It can then 
be shown [16, 19| that if a resummation of the form of Ref. holds {i.e., of the form Eq. ( 4.61 ) 
with all Qij = when i > 0), the MS anomalous dimension 7"^^ has, to all orders in a^, at most 
a single ln(l/A^'^). In order to get from the physical to the MS scheme one has to perform a 
suitable combination of the well-known unresummed scheme change which takes to MS from 
the so-called DIS scheme pOf, and a resummed scheme change. This latter resummed scheme 
change is in fact entirely determined by the requirement that the MS anomalous dimension be 
linear in ln(l/A^"). 

To see this, assume that the resummed anomalous dimension in the physical scheme is given 

by 



7 



,AP,i 



and perform a resummed scheme change Eqs. ( |6.12| , |6TT3D , with 



(6.15) 

for i > 0. After such a 



scheme change 

/AP,res 



7'^^'^^^(Ar,a.(;.^)) 



n j=i 



n j=i 



(9oj + zoj) ai{^?/n) + In — ^ ^a^/i^ 



n 



(6.16) 



Hence, imposing that 7' 



AP.res 



be linear in ln(l/A^") fixes 

00 

Z^^^ = -Y.g,,ai{k^/n). 



(6.17) 



Performing the scheme change Eq. ( |6.17| ) one gets 



7 



/AP.res / 



(iV,«s(/i')) 

(Ar,l,a,(g 



^ /-^"rfn rQ' dk^ ,,,2 , , 



(6.18) 
(6.19) 
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The resummed MS anomalous dimension and coefficient function can now be obtained by trans- 
forming those given in Eqs. ( |6 .181 , 1091) with an unresummed scheme change obtained by using 
of Ref. but omitting all contributions to Z-^ which are proportional to ln™'(l/iV'^) with 
m > 1, so that the requirement on discussed after Eq. ( |6.14| ) is satisfied. 

In practice, however, it is easier to observe that the requirement that the MS anomalous 
dimension be free of higher order powers of In simply means that the MS resummed anomalous 
dimension coincides with the unresummed one. This further implies that the MS resummed 
coefficient function is found combining the resummed result Eq. (|6.19| ) with the standard MS 
unresummed coefficient function, and subtracting double counting. At next-to-leading order, 
next-to-leading log we get 



Cjj^[N,l,as{Q^)) = exp 



exp 



1 n JQir 



f^^" dn f^- dk' f .,2 / N 2/,2 / N 

/ — / -JY [%i(^s[k In) +go2a,{k jn) 

Jl n JQ^n K ^ 

fO(af+2ln^Ar) +0(«2) 

p (goiasik'^/n) + go2al{k'^/n)^ 

+ as{Q'){c^iN) + ^\n'^^+0{a':^Hn'N)+0{al) (6.20) 

where Ci{N) is the standard next-to- leading order contribution to the coefficient function. In the 
first line on the r.h.s. we have given the result which is obtained by performing the two scheme 
changes Z^Z^, which correspond respectively to the first and second factor in square brackets, 
and in the second line the result obtained by combining Eq. (|6.19|) with the unresummed result 
and explicitly performing the double-counting subtraction; of course the two results agree up to 
subleading corrections. 

Resummed results in the MS scheme can be analogously obtained starting from the expression 
Eq. ( |4.63| ) of the resummed physical anomalous dimension in terms of a resummed splitting 
function. Proceeding in a similar way, with obvious modifications, we get 



C^(N,l,as{Q' 



X 



,N~1 




-p- [goias{k ) +go2a,{k 



+ 0{a':+Hn''N)+0{al), (6.21) 



where now goi are contributions to the coefficient of the 0(ln(l — x)") term in the physical 
splitting function, and Ci(x) is the x-space next-to-leading order contribution to the coefficient 
function. Note that the double-counting subtraction in Eq. ( |6.20|) is not just the Mellin transform 
of the double counting subtraction of Eq. ( |6.21| ): this is only true at the leading log level. Again, 
the two forms Eq. (|6.20|) and ( |6.21|) of the resummed MS coefficient function coincide at the 
next-to-leading order and next-to-leading log level, but differ by terms which are at least of 
0{al) and then either next-to-next-to- leading log or suppressed by powers of A^. 

The resummed result Eq. ( |6.21| ) can be finally brought in the form of Ref. by recalling 
that §02 is determined in terms of the physical anomalous dimension, and separating off the 
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contributions which originate from the anomalous dimension and coefficient function, according 
to Eq. ( p36|) : 



(6.22) 



/•I , x^~^ - 1 

exp a / ax — ; 

JO I — X 



)2 A;2 



+asiQ-) I X 



N-l 



1 — X 



h B^as{Q\l - xT) 



Note that the double counting subtraction is not the same as in Eq. ( |6.21| ), because the 5-term 
in Eq. ( |6.22|) is not free of 0{lnj^) terms. 

The physical and MS schemes which we discussed explicitly are extreme choices in the treat- 
ment of the resummed terms: in the former case, these are entirely included in the anomalous 
dimension, and in the latter case they are entirely included in the coefficient function. Of course, 
a variety of intermediate choices is possible. These resummed scheme choices can then be freely 
combined with the usual choice of unresummed factorization scheme. For instance, one can 
combine an MS-like resummation (where the resummation is included in the coefficient func- 
tion) with a physical-scheme anomalous dimension. In such case, the anomalous dimension is 
given, up to 0{N^) terms, by Eq. ( |6.1(j| ). Clearly, all such scheme choices can be obtained by a 
combination of the resummed and unresummed scheme changes Eqs. (|6.9|,|6.1lD. 



7 Summary 

In this paper, we have discussed the resummation of large logarithmic contributions to pertur- 
bative coefficients that arise as a left-over of the cancellation of infrared singularities near the 
boundary of the phase space. Our approach is based on identifying the relevant dimensionful 
scale for the large logs, showing that it is determined essentially by kinematic considerations, 
and using the renormalization group to resum it. 

The main result of this paper is a resummation formula which allows the computation of 
the physical anomalous dimension {i.e., the scale dependence of physical observables) to any 
logarithmic accuracy, starting from a fixed-order calculation. This result is less predictive than 
what one gets by extending the structure of familiar next-to-leading log resummation |Q, |^ 
to high orders: specifically, we find that an order M^hii perturbative calculation is required for 
resummation to next''~^-to-leading log accuracy. In particular, the standard next-to-leading log 
resummation turns out to be fully determined by the fixed next-to-leading order result only 
because of the vanishing of a coefficient which, in our approach, appears to be accidental. 

The advantage of this approach, however, is that it does not rely on factorization properties 
of the physical cross section. Indeed, our proof does not require any detailed consideration of the 
individual Feynman diagrams which contribute in the soft limit, and thus in particular it does 
not require the classification of the different forms of radiation which contribute in the soft limit 
(soft radiation, 'jet' radiation coUinear to incoming or outgoing partons etc.) which is the main 
complication in conventional next-to-leading 0, |^ or all-order approaches to resummation. 
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An important feature of the present approach is that resummed results are provided directly 
for a physical observable, the physical anomalous dimension, rather than for partonic quantities. 
This makes it easier to study issues related to the choice of an appropriate resummation variable. 
In particular, we have shown that a leading (or next-to-leading, etc.) ln(l — x) expression of 
the cross section is necessarily ill-defined, and traced the origin of this problem, first pointed 
out in Ref. 0, to subleading ln(l — x) terms which appear in the inverse Mellin transform of 
a leading InA^ expression, which is instead well-defined. Also, the discussion of issues related 
to the choice of factorization scheme is considerably streamlined: our results are manifestly 
scheme-independent, and results within specific choices of scheme can be readily obtained as 
special cases. 

Because of its generality, the approach to soft resummation presented in this paper lends 
itself naturally to a variety of extensions and generalizations. In particular, it should be easy 
to generalize it to other classes of soft resummation, such as jet or prompt-photon production, 
and, perhaps more interestingly, to the resummation of transverse momentum distributions, 
which can be unified Q with the resummation discussed here by considering Mellin and Fourier 
transforms as special cases of a complex integral transform. 

It is interesting to ask whether the formalism presented in this paper might also lead to 
further insight on resummation at small x. 
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A N-body phase space 

The n-body phase space can be expressed in terms of the n — 1-body and the two-body phase 
spaces To see this, start with the standard definition of the phase space for a generic process 
with incoming momentum P and n bodies in the final state with outgoing momenta ki,. . . ,kn 
and masses kf = mf: 

dUP; h,..., = • • • (^$^(2vr)'^5('^)(P - fc, - ... - k^), (A.1) 

The momentum-conservation deha function can be rewritten as 

Now note that 
where 

M^i = {Pn-i? - \Pn-i?; (mi + m2 + . . . + m„_i)' < M^^ < (v^ - m„)' . (A.4) 
We obtain immediately 

j-s 

d(j)n{P] ki,..., kn) = ~(2^ #2(^; Pn-l, ^n) d(j)n-l{Pn-i; ^1, • • • , kn-l), (A.5) 

which is the desired result. Using this result recursively, the n body phase space can be entirely 
expressed in terms of two-body phase space integrals. 

For completeness, we also give the expression of the two-body phase space in the center-of- 
mass frame. We have 

dMP.^uP.,h) = (2.Plp;o i2fr-'^2k? (^")'^^'^(^-+^ - ^' - 

= ^-^-p^HP'L-Prk°). (A,6) 

In the center-of-mass frame, Pj+i = 0, P^^^ = Mj+i and |Pj| = \ki\ = k^, since kf = 0. We have 
therefore 

Hence 

\k-\ -k^- - -Ml±^(i_J^\ (A 8) 

(A.9) 



and 
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